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ABSTRACT 

We report the discovery of a 200 mHz quasi-periodic oscillation (QPO) in the X-ray emission from a 
bright ultra- luminous X-ray source (ULX) Holmberg IX X-l using a long XMM-Newton observation. 
The QPO has a centroid at vqpo — 202.5^|'g mHz, a coherence Q = vqpo / ^ v fwhm ~ 9.3 and an 
amplitude (rms) of 6% in the 0.2 — 10 keV band. This is only the second detection of a QPO from 
an ULX, after M 82 X-l, and provides strong evidence against beaming. The power spectrum is well 
fitted by a power law with an index of w 0.7. The total integrated power (rms) is ss 9.4% in the 
0.001 — 1 Hz range. The X-ray spectrum shows clear evidence for a soft X-ray excess component that 
is well described by a multicolor disk blackbody (fcTj„ ~ 0.3 keV) and a high energy curvature that 
can be modeled either by a cut-off power law (T ~ 1; E cuto ff = 9 keV) or as a strongly Comptonizcd 
continuum in an optically thick (r « 7.3) and cool (kT e « 3 keV) plasma. Both the presence of the 
QPO and the shape of the X-ray spectrum strongly suggest that the ULX is not in the high/soft 
or thermally dominated state. A truncated disk and inner optically thick corona may explain the 
observed X-ray spectrum and the presence of the QPO. 

Subject headings: accretion, accretion disks — stars: individual (Holmberg IX X-l; M81 X-9) - 
X-rays: stars 



1. INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are extra- 
nuclear point X-ray sources with luminosities exceed- 
ing the Eddington limit for a 10 M Q black hole (BH). 
The most popular model to explain the high luminosi- 
ties of ULX is the "intermediate mass BH" with mass 
M BH ^ 10 2 -10 4 M Q (e.g., Colbert & Mushotzky 1999). 
Other popular models include X-ray binaries (XRBs) 
with anisotropic emission (King et al. 2001), beamed 
XRBs with relativistic jets directly pointing towards us i. 
e., scaled down versions of blazars (Mirabel & Rodriguez 
1999), and XRBs with super-Eddington accretion rates 
(Begelman 2002). 

There is now extensive study of X-ray spectra of ULXs 
with XMM-Newton & Chandra. ULXs show a vari- 
ety of spectral shapes: (i) simple power law similar to 
the low/hard state of BH XRBs (Winter, Mushotzky & 
Reynolds 2005), (H) disk blackbody (kT ~ 0.1-0.4 keV) 
plus power law similar to the high/soft state of BH XRBs 
(see e.g., Feng & Kaaret 2005), (Hi) strongly curved spec- 
tra at high energies (Stobbart, Roberts & Warwick 2004; 
Feng & Kaaret 2005; Agrawal & Misra 2006; Dewangan, 
Griffiths & Rao 2006a; Stobbart, Roberts & Wilms 2006) 
without a corresponding spectral state in BH XRBs. 
These spectra are physically well described by a cool disk 
and thermal Comptonization in an optical thick corona. 
The cool disk plus power-law spectra of ULXs may not 
always correspond to the high/soft state of BH XRBs (see 
Roberts et al. 2005; Dewangan et al. 2006a). The de- 
tailed variability properties of ULXs are yet to be known. 
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It is important to study both the spectral and temporal 
characteristics of ULXs to understand their nature. BH 
XRBs exhibit characteristic power density spectra (PDS) 
depending on their X-ray spectral state (McClintock & 
Remillard 2003). Many BH XRBs show quasi-periodic 
oscillations (QPOs) and breaks in their PDS that rep- 
resent characteristic timescales close to the BH. These 
characteristic frequencies scale with the mass of the BH. 
Therefore, determination of the shape of the PDS and 
the detection of QPOs is crucial to understand the na- 
ture of ULXs. To date, there is only one ULX, M 82 X-l, 
that is known to show QPOs (Strohmayer & Mushotzky 
2003; Fiorito & Titarchuk 2004; Dewangan, Titarchuk & 
Griffiths 2006b; Mucciarelli et al. 2006). 

In this Letter, we report the discovery of a 200 mHz 
QPO from the bright ULX Holmberg IX X-l (hereafter 
Ho IX X-l; also known as M81 X-9) based on a long 
XMM-Newton observation. The ULX is located in the 
dwarf irregular galaxy Holmberg IX, a companion to 
M 81. The position of the ULX is ~ 2' away from 
the optical center of the host galaxy. Ho IX X-l is a 
bright X-ray source (Lx ^ 10 40 erg s -1 ) that is variable 
on time scales of weeks and months (La Parola et al. 
2001). Miller, Fabian & Miller (2004) presented X-ray 
spectra of Ho IX X-l based on two snapshot XMM- 
Newton observations. The spectra were modeled by a 
cool (kT ~ 0.17 — 0.29) disk emission and simple power- 
law components. Miller et al. (2004) suggested a BH 
mass of 10 3 M Q by scaling the disk temperatures mea- 
sured in Ho IX X-l to those usually observed in BH 
XRBs in their high state. Here we present power and 
energy spectral analysis and show evidence for a QPO 
and high energy spectral curvature. 

2. OBSERVATION & DATA REDUCTION 

XMM-Newton observed Holmberg IX on 2004 Septem- 
ber 26 for 119 ks . The EPIC PN and MOS cameras were 
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operated in the large window and the full frame modes, 
respectively, using the thin filter. We used SAS version 
6.5.0 and the most recent calibration database to pro- 
cess and filter the event data. Examination of the back- 
ground rate above 10 keV showed that the observation is 
completely swamped by the particle background after an 
elapsed time of 105 ks and this latter period was there- 
fore excluded from the rest of the analysis. For temporal 
analysis, we used all the PN and MOS events with pat- 
terns < 12 and a continuous exposure of 103.7 ks during 
which both the PN and MOS cameras operated simul- 
taneously. For spectral analysis, we chose a count rate 
cut-off criterion to exclude the high particle background 
and used event patterns — 4 (singles and doubles) and 
FLAG=0. We extracted events from the PN and MOS 
cameras using a 45" circular region centered at the peak 
position of the ULX {a = 09 h 57 m 54 s , S = 69°03'46"; 
La Parola et al. 2001). We also extracted appropriate 
background events within nearby circular regions free of 
sources. 

3. THE POWER SPECTRUM 

For temporal analysis, we combined the PN and MOS 
data to increase the signal-to-noise ratio. We calculated 
a power density spectrum (PDS) using the background 
corrected PN+MOS light curves sampled at 0.5 s. Fig- 
uredshows the 0.2 - 10 keV PDS of Ho IX X-l rebinned 
by a factor of 1024 yielding a frequency resolution of 
7.8 mHz. The PDS continuum falls off approximately 
following a power law up to ~ 60 mHz where the white 
noise arising from the Poisson errors starts to dominate 
the continuum. The white noise continuum above 0.5 Hz 
is not shown in Fig. ^ There is a prominent QPO with 
its peak frequency near 200 mHz. 

To fit the PDS, we used a simple power law for the con- 
tinuum, a Lorentzian to describe the QPO and a constant 
to account for the white noise. This model resulted in 
an acceptable fit, providing a minimum \ 2 = 28.4 for 58 
degrees of freedom (dof ) . In Fig. ^ we show the best- 
fit model to the PDS as a thick line. The errors on the 
best-fit PDS model parameters, quoted below, are at a 
la level. The best-fit model resulted in a QPO centroid 
frequency vqpo = 202.6 ± 3.8 mHz, a width vfwhm = 
21.8^° 7 5 mHz and an amplitude A QPO = (3.7 ± 0.9) x 
10~ 3 (rms/mean) 2 . The best-fit power-law index and 
the amplitude are V = 0.73 ± 0.02 and A PL = (7.6 ± 
0.7) x 10~ 3 (rms/mean) 2 / Hz at 0.1 Hz, respectively. 
The best-fit constant is 0.826 ±0.005 (rms/mean) 2 / Hz. 
The total integrated power (0.001 - 1 Hz) and the QPO 
power expressed as rms/mean are 9.4% and 6%, re- 
spectively. To estimate the statistical significance of 
the detection of the QPO, we also fitted the PDS with 
a simple power law and a constant. This model pro- 
vided a minimum x 2 = 41.5 for 63 dof. Thus the addi- 
tion of the Lorentzian for the QPO was an improvement 
(A% 2 = —13.1 for three additional parameters) at a sig- 
nificance level of 99.5% based on the maximum likelihood 
ratio test. 

4. THE ENERGY SPECTRUM 

The PN and MOS spectra were grouped to a minimum 
of 50 and 20 counts per spectral channel, respectively, 
and analyzed with XSPEC 11.3. The errors on the best- 
fit spectral parameters are quoted at a 90% confidence 
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Fig. 1. — Power density spectrum of Holmberg IX X-l derived 
from the EPIC PN and MOS data in the 0.2-10 keV band. The 
white noise level expected (~ 0.8) from the Poisson errors has 
not been subtracted. The frequency resolution is 7.8 mHz. The 
best-fitting model comprising a simple power-law for the PDS con- 
tinuum, a constant for the white noise and a Lorentzian for the 
QPO is shown as a thick line. 
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FlG. 2.— Ratios of EPIC PN data and the best-fit absorbed PL, 
MCD+PL and MCD+cutoff PL models. The soft X-ray excess and 
the high energy curvature are clearly seen in the top panel. The 
data have been rebinned for the purpose of display only. 



level. First we fitted a simple absorbed power law (PL) 
model to the PN and MOS spectra of Ho IX X-l. We fit 
the MOS1 and MOS2 data jointly with an overall nor- 
malization constant to account for possible differences in 
the source extraction areas or calibration uncertainties. 
We used the 0.2 — 12 keV band in all the fits. The simple 
power law model resulted in minimum \ 2 — 1850.3 for 
997 dof and 2365.1 for 1906 dof for the PN and MOS 
data, respectively, thus providing statistically unaccept- 
able fits to both the PN and MOS data. We have plotted 
the ratio of the PN data and the best-fit PL model in Fig- 
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Fig. 3.— The EPIC PN data and the best-fit absorbed MCD 
plus thermal Comptonization model. Also shown is the ratio of the 
PN data and the best-fit model. 



TABLE 1 

Best-fit spectral model parameters for Holmberg IX X-l. 



Parameter 


Model A a 


Model B a 




PN 


MOS 


PN 
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N H h 


1 85+ ' 12 


i 09+0.10 

- L, ° -0.0.07 


-, 7S +0.12 
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1 26+ - 2 
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'--'edge 


66+ 02 
u - DD -0.02 
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u - DD -0.02 
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u - 1D -0.04 




16+ 003 
u - 1D -0.02 
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u - lo -0.07 
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n 00+O.O2 
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29+ 04 
u,ZM -0.03 


A M CD d 


on 1+4.2 
oa ' -3.0 
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i2.o±f;| 


r 


1 02 +0 10 
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'• -Q.5 


9 9+ 04 


E Ka c 


f. 99 +0.03 
°' z -0.05 
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0.01 (f) 




0.01(f) 






1 4+0-8 

^-o.s 




1 4+0-8 




EW Ka c 


0.027 




0.026 




fx 1 


6.0 


6.1 


6.0 


6.1 


Lx ! 


8.3 


8.4 


8.3 


8.4 


X 2 - 

A mm 

dof 


983.3 


1965.0 


982.3 


1967.4 


988 


1903 
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a Model A: 2Edge(MCD+cutoff PL+Gaussian); Model B: 
2Edge(MCD+compTT+Gaussian). For the MOS data, the edges 
and the Gaussian line were excluded from the above models. 

b In units of 10 21 cm -2 . 

c In units of keV. 

d The MCD normalization A MCD = (R in / km) / (D / 10 kpc), 
where Ri„ is the inner radius and D is the dis- 
tance. Apl is the power-law normalization in units of 
10 -4 photons cm -2 s _1 keV -1 at 1 keV. 

c In units of 10 — 6 photons cm -2 s~ x . 

'Observed flux and luminosity in units of 10 -12 ergs cm -2 s — 1 
and 10 39 erg s — 1 , respectively, and in the band of 0.2 — 10 keV. 



ure[5](top panel). This plot clearly show a broad hump in 
the 3 — 8 keV band and a soft X-ray hump or excess emis- 
sion below 1.5 keV. The shape of the X-ray spectrum of 
Ho IX X-l is very similar to that observed for NGC 1313 
X-l (Dewangan et al. 2006a). Similar ULX spectra have 
also been reported by Feng & Kaaret (2005) and Sto- 
bbart et al. (2006). These spectra are well described 
by a steep power law at low energies and a hot accre- 
tion disk emission at high energies (Stoobart et al. 2006; 
Feng & Kaaret 2005) or by a cool disk component and a 
cutoff power-law (Dewangan et al. 2006a). Addition of 
a hot multicolor disk blackbody (MCD) component im- 
proved the fit to the PN data significantly (x 2 — 1149.5 
for 995 dof). Significant improvement is found for the 
MOS data also. Fig. ^(middle panel) shows the ratio of 
the PN data and the best-fit MCD+PL model. The high 
energy curvature is well described by the hot MCD com- 
ponent but the soft excess component is not well modeled 
by the steep power-law component. Replacing the sim- 
ple power-law with a cut-off power law improves the fit 
(PN: x 2 = 1028.6 for 994 dof and MOS: x 2 = 1965.0 
for 1903 dof) and the MCD component results in cool 
disk (kTi n ~ 0.25 keV). Examination of the ratio of the 
PN data and the best-fit MCD+cutof f PL model, shown 
in Fig. [21 {bottom panel), reveals several weak features 
e.g., possible edges at ~ 0.7 and ~ 8 keV and an emis- 
sion feature at ~ 6 keV. Addition of edges at ~ 8 keV 
and ~ 0.7 keV to the MCD+cutoff PL model improves 
the fits to the PN data by Ax 2 = -8.0 and -26.2, re- 
spectively, for two additional parameters. Addition of 
a narrow Gaussian line at ~ 6.2 keV improves the PN 
fit marginally (A% 2 = —11.0 for two additional parame- 
ters). The best-fit parameters are listed in Table The 
MOS data do not show evidence for any discrete feature, 
addition of edges or lines does not improve the fit. 

The high energy curvature or cut-off around ~ 7 keV 
can physically be described as arising from thermal 
Comptonization in an optically thick corona (Goad et al. 
2006; Agrawal & Misra 2006; Stobbart et al. 2006). To 
test such a scenario for Ho IX X-l, we used the XSPEC 
model compTT (Titarchuk 1994) that describes thermal 
Comptonization of soft photons in a hot plasma. The pa- 
rameters of this model are the soft photon temperature 
(kTs), electron plasma temperature (kT e ), plasma opti- 
cal depth (r) and the normalization. We also used the 
MCD model to describe the soft excess component. The 
MCD+compTT model provided minimum % 2 values similar 
to that obtained for the MCD+cutoff PL model. Addi- 
tions of two absorption edges at ~ 8 keV and ~ 0.7 keV 
and a narrow Gaussian line at ~ 6.2 keV improved the 
fit to the PN data by A X 2 = -8.8, -30.1 and -8.2, 
respectively. The best-fit parameters are listed in Ta- 
ble n The curvature or cut-off is well modeled by ther- 
mal Comptonization in an optically thick (r = 7.3^ 'g 

and relatively cool (kT e = 3.01q4 keV) electron plasma. 
Using the above model, we derive an unabsorbed flux 
of 1.06 x 10~ n ergs cm -2 s _1 for Ho IX X-l in the 
0.1 — 100 keV band. This translates into a bolometric 
luminosity of ~ 1.5 x 10 40 ergs s _1 using a distance of 
3.4 Mpc (Georgiev et al. 1991; Hill et al. 1993). 

5. DISCUSSION 

Based on a long XMM-Newton observation, we have 
discovered a QPO with centroid frequency vqpo — 
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202. 5lg 8 mH z in the power spectrum of Ho IX X-l. The 
QPO has an amplitude (rms) of 6% and a coherence 
Q « 9.3. The underlying PDS is well described by a 
power law with an index of « 0.7. This is only the second 
detection of a QPO from any ULX. M82 X-l was the first 
ULX to show clear evidence for a QPO (Strohmayer & 
Mushotzky 2003). M82 X-l shows a variable QPO with 
the centroid frequency in the range of « 50 — 170 mHz 
(Muccarelli et al. 2006). The 200 mHz QPO from 
Ho IX X-l appears to be similar to the 54 mHz QPO (rms 
amplitude « 8.5%; Q « 5; Strohmayer & Mushotzky 
2003), seen in the first XMM-Newton observation of M82 
X-l, in terms of the amplitude and coherence. The PDS 
continua of Ho IX X-l and M 82 X-l (in the first XMM- 
Newton observation) are also simple power laws without 
a clear evidence for a break. 

BH XRBs show flat-topped, band-limited noise in their 
low state (LS). The flat-topped noise weakens in the in- 
termediate/very high states (IS/VHS) and is absent in 
the high state (HS). Strong low frequency QPOs are ob- 
served in the IS /VHS and occasionally in the LS but are 
usually not present in the HS (see e.g., McClintock & 
Remillard 2003). The IS/VHS are also known as the soft 
or hard intermediate states (SIMS or HIMS), depending 
on whether the spectrum is soft or hard, in the classifi- 
cation scheme of Belloni (2005). In the HIMS, BH XRBs 
show strong 0.1 - 15 Hz QPO with high Q-factor (7-12) 
and large amplitude (3 — 16% rms). These are called type 
C QPOs (Casella et al. 2005). The QPO in Ho IX X-l 
(as well as in M82 X-l) have properties similar to the 
type C QPOs and suggest that the variability properties 
of these ULXs are likely to be similar to that of BH XRBs 
in their HIMS. 

X-ray spectrum of Ho IX X-l is curved at high en- 
ergy (E ^ 6 keV) and shows a soft X-ray excess emission 
at low energies ;$ 1 keV. The PN spectrum is well de- 
scribed by an MCD (kT ~ 0.25 keV) and a cut-off power 
law (r ~ 1; E cuto ff ~ 9 keV) or a compTT component 
(kT s ~ 0.5 keV; kT e ~ 3 keV and r ~ 7). While this 
form of unusual X-ray spectrum is common among bright 
ULXs (Stobbart et al. 2006), there is no corresponding 
spectral states of BH XRBs. Both the presence of the 
QPO and the high energy curvature in the X-ray emis- 
sion from Ho IX X-l are inconsistent with the ULX being 
in the high/soft state and the soft X-ray excess emission 
cannot be interpreted as the optically thick emission from 
a disk extending down to the last stable orbit. Thus the 
temperature of the soft excess component cannot be used 
to estimate the BH mass for Ho IX X-l. 

In BH XRBs, QPOs are observed when both the disk 
and power-law components contribute significantly to the 
total X-ray emission (e.g., Sobczak et al. 2000). When 
strong QPOs are seen in the LS, IS and VHS, the tem- 
perature of the disk component is cooler than that in the 
HS. In these cases, the inner disk is thought to be trun- 
cated (Done & Kubota 2005). The presence of the strong 
QPO in the power spectrum of Ho IX X-l, then, implies 
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that its accretion disk must be truncated at some inner 
radius, larger than the radius of the innermost stable or- 
bit. Done & Kubota (2005) have investigated accretion 
disk structure of BH XRBs in their VHS characterized by 
strongly Comptonized spectra and cooler disk tempera- 
tures compared to that in the HS. They have proposed 
a model in which the corona and the disk are parts of 
the same accretion flow. At high accretion rates, the in- 
ner corona is energetically-coupled to an outer disk that 
is truncated at some inner radius, thus resulting in a 
cooler disk and optically thick inner corona. ULXs with 
cool soft X-ray emission and high energy curvature such 
as Ho IX X-l may be the extreme cases of VHS of BH 
XRBs. In Ho IX X-l, the disk truncation radius is likely 
to be larger and the inner corona is optically thicker com- 
pared to that in the VHS of BH XRBs. Thus Ho IX X-l 
may be accreting at a rate (relative to the Eddington 
rate) larger than that in the VHS of BH XRBs. Interest- 
ingly, some BH XRBs e.g., GRS 1915+105 show spectral 
curvature around 30 keV in their HIMS (Vadawale et al. 
2001; Kubota & Done 2004). Due to the presence of a 
strong QPO and spectral curvature at ~ 3 keV in ULXs, 
Ho IX X-l and M 82 X-l, similar to the type-C QPO 
and high energy curvature at ~ 30 keV in the HIMS of 
some BH XRBs, it is quite conceivable that ULXs have 
their spectral curvature at lower energies in their HIMS. 

If the association of the observed QPO in ULXs to 
type-C QPO of BH XRBs is correct, it gives a useful 
handle to constrain the mass of the BH in Ho IX X- 
1. In the HIMS states, the accretion rate is typically 
> 30% of the Eddington luminsoity. Considering the 
fact that QPOs are associated with accretion discs (and 
hence ruling out the possibility of beamed jet emission for 
the observed X-ray emission), this puts a strong upper 
limit of 400 M Q for the mass of the BH in Ho IX X- 
1. Further, the spectral curvatures in GRS 1915+105 
are preferentially seen just before radio flares and the 
QPO frequencies during such states are in the range of 
0.7 — 3 Hz (Vadawale et al. 2003). If we equate this 
frequecy range to the observed QPO in Ho IX X-l, we get 
an estimate for the BH mass in Ho IX X-l as 50—200 M Q . 

Associating the observed QPO in Ho IX X-l to type C 
QPO in BH XRBs has another interesting consequence. 
It is observed that transition from HIMS to soft state 
(or SIMS) gives rise to fast relativistic jet in galactic 
sources. Similarly, such ultra-luminous jets should be 
seen in ULXs, when they make a transition from HIMS 
to SIMS. Detecting such relativistic jets in ULXs will put 
the association with the BH XRBs on a firm footing. 
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